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I. INTRODUCTION 



The medium induced energy loss of a hard (quark or gluon) jet traversing matter - hot or cold - has recently been 
the subject of intensive interest pj-p^. The radiative energy loss has been found to be independent of the jet energy, 
for large energies, and growing as L^where L is the extent of the medium |^,^. The order of magnitude of this effect 
in hot matter may be expected to be large enough - compared to the case of cold nuclear matter - to lead to an 
observable and remarkable signal of the production of deconfined matter. 

The natural observable to measure this energy loss would then be the transverse momentum spectrum of hard jets 
produced in heavy-ion collisions. Jet quenching which has been discussed recently in |l^Jl^ is the manifestation of 
energy loss as seen in the suppression and change of shape of the jet spectrum compared with hadronic data. One may 
also think of measuring single particle inclusive spectra. Comparisons of high pj- particle spectra in p + p,p + A and 
A + B collisions at SPS energies have already been performed and discussed [Q. The results seem to be difficult to 
interpret in view of the obvious model dependence of the theoretical prediction. This model dependence is associated 
with the relatively low energy and p± ranges which are currently accessible fl^ . The situation should be much more 
favorable at RHIC. Jet production will be intensively studied at LHC. 

In Section 2 of the present paper we study the angular distribution of radiated gluons which are believed to be 
the main source of energy loss. This allows for quantitative predictions of the energy lost outside the cone defining 
the jet. This problem has already been considered in |p|-[rT[|, where the expression of the "characteristic" angle for 
gluon emission has been worked out. We work along the same lines but we perform a more complete calculation of 
the angular distribution. We concentrate on the realistic case of a hard jet produced in the medium. We make the 
implicit assumption that the lifetime of the hot medium is large enough to allow for a large number of scatterings of 
the jet as it traverses the medium. 

In Section 3 we calculate the integrated energy loss outside an angular cone with fixed opening angle, AE(9cono), 
and derive the expression for the ratio R{9conc) = Ai?(6'cono)/A£^, where AE is the completely integrated loss. For 
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simplicity of expression we take a convention where Ai? is a positive quantity. The surprising result is that i?(6'conc) 
is a universal function of the variable O^^^^qL^ , where q is the transport coefficient characteristic of the medium. 

In Section 4 we give estimates and in particular we show that the energy loss is more coUimated in the case of a 
hot QCD medium as compared to nuclear matter. As a consequence the characteristic angle for gluon emission in a 
hot medium is quite small of the order of 10° for a temperature of 250 MeV. It is, however, possible to choose large 
enough angles of the order of 30° and still have an appreciable energy loss. We also make a comparison with the 
corresponding energy loss outside a cone in the vacuum. All the calculations are done on the partonic level with no 
attempt to include fragmentation of partons into hadrons. Fragmentation effects, as well as effects due to plasma 
expansion iQ, could be important and should be studied. 

Technical details are summarized in Appendices A - C. 



II. MOMENTUM SPECTRUM OF RADIATED GLUONS 



In this section we discuss the induced momentum spectrum of soft gluon emission {x 0) from a fast quark jet. We 
assume that the quark is produced inside the medium by a hard scattering at time t = 0. From the production point 
it propagates over a length L of QCD matter, carrying out many scatterings with the medium. We follow closely the 
derivation and notation given in BDMPS |,| (denoted by I and II) and EDMS @ which we denote by III. 

In the case under consideration the spectrum per unit length z of the medium consists of two terms: one correspond- 
ing to an "on-shell" quark, as if that quark were entering the medium, and one corresponding to a hard production 
vertex as described in III (cf. Eq.(31)): 

Ludl LudI 
dLjdzd'^lJ_ dujdzd'^U_ 

uj denotes the soft gluon energy and U_ the scaled transverse gluon momentum U_ — k/^ with /i the appropriate scale 
of the QCD potential. 

As follows from Eq.(31) in III the induced gluon spectrum is 

(2) 

k 

where fh differs from / in that the gluon emission time has been evaluated at the time of the hard interaction, t ~ 0. 
After we have converted to impact parameter space the explicit relationship between and / will be given. In 
the soft CO limit the contributions from the diagrams (Fig. la-c) are included. In comparison with the diagrams shown 
in Fig. 4 of III we note that the soft limit leads to a vanishing of the sum of diagrams Fig. 4d-4j. 
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FIG. 1. The diagrams representing the contributions to the spectrum (g). 

The quantities in (|^), explained in III, have the following meaning: agCp /t^'^ is the coupling of a gluon to a quark. 
The factor -^^fiU. + Q,t2 — ti)padti gives the number of scatterers in the medium, padti, times the gluon emission 
amplitude at ti, evolved to t2. The momentum labels are shown in Fig. 1. We describe the scattering in terms 
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of the normalized cross section V{Q) = -^-^^ which depends on the scaled momentum Q = q/ ^i. We note that 
the final state gluon carries a transverse momentum k. The factor corresponding to the amplitude /q in III which 



sums the gluon emissions (Fig. la-c) is / (PQ V{Q) 2 



U+Q 



U 

IF 



In order to eliminate the medium-independent 

^ . The parameter 



_(F+qF 

factorisation contribution one has to perform a subtraction of the value of the integrals at k = 

K 

, 2Cf Xfi^ 



(3) 



depends on medium properties, in particular on the quark's mean free path A — 1/ pi, as well as on the gluon energy. 

In dealing with (^) in contrast to the [/-integrated spectrum discussed in III we now have to take into account the 
possibility that the emitted gluon may rescatter in the medium. For example the contribution shown in Fig. la may 
now have final state interactions as illustrated in Fig. 2. 




FIG. 2. Final state interactions contributing to Fig. la; (a) real, (b) and (c) virtual interactions. 

This is taken care of by the final state interaction amplitude Tfsi in @j which is most conveniently expressed in 
impact parameter space. One goes to impact parameter space by taking the Fourier transform, for example for the 
ampHtude /(f/, i), 

/(C/,<) = / ^,e^^-^f{B,t). (4) 

As in III we rescale the time variable t — ^^Ar and define 

To = N,L/2Cf\. (5) 
Following our analysis of parton px-broadening in II we write the final state interaction amplitude as 

1 



Tfsi{B^ ~ B2,t) = exp 



^{B^'B^Ytv 



(6) 



In II the distribution in transverse momentum and the corresponding distribution in impact parameter were given 
for a single parton going through a medium. Here we have a quark-gluon system passing through the medium. 
However, interactions of the quark with the medium cancel between real and virtual interactions leaving only the 
interactions of the gluon with the medium as in II. These interactions can be real interactions, as shown in Fig. 
2a where both the gluon in the amplitude at impact parameter B_^ and in the complex conjugate amplitude at Bj2 
interact with the medium, or they can be virtual interactions as shown in Figs. 2b and 2c. The real interactions give 
a factor V{B_-y ~ B.2) while the virtual interactions give the absorption of the gluon wave and a factor of — V^(0) = — 1. 
V{B) is the Fourier transform of the potential V{Q). These terms combine together giving at small B_i ~ B_2 

l^ViB,-B,)^^iB,-B2)^v. (7) 

V has a logarithmic dependence on B^. In the following we neglect the variation of the logarithm, an approximation 
which is good except when large transverse momentum rescatterings occur (see II) . The factor in (^ then exponenti- 
ates, as in II, to give (H), where the 5 in the exponent (||) differs from the j in Eq.(pl|) of II because of the difference 
of T from t used there. Using (||) and (^) in (||) wc obtain the spectrum in terms of impact parameter integrals. 
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LudI 



dtodzcPU 
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'iTTiB_2 



V{B,-B2)-V{B,) 



K=0 



(8) 



We have used fh {B_, r) 



vB- 



■f{B_, t) which foUows from the fact that fh and / have the same time evolution while 



the initial conditions satisfy fh{B_, 0) — 0) as can be seen by comparing Eqs.(16) and (17b) of III in the soft 

gluon limit. 

When integrating (^) with respect to the gluon momentum U_ it is straightforward to reproduce the energy spectrum 
given e.g. in Eq.(34) of III. Next we recall that the amplitude f{B_, r) is given by (cf. Eq.(38) in III) 



fiB,r) 



COS^ WqT 



Bexp —-mujoB_ tanwoT , 



(9) 



with 



'1/2k and ljo = \ 



(10) 



together with k given in (||). In this approximation (leading logarithm and small a;) after substituting (^ into (^ we 
find a convenient expression for the momentum spectrum per unit length 



-V Re 



diodzd'^Li L 
with the time dependent variables 



Jo COS^WoT [Jq 7-00 J 



K=0 



(11) 



i vtsinujoT 
a — —mujQ tantJoT = — - 



2 luq 



/3'-^(To-r). 



(12) 



After using (0) in (||) the B-space integrals become 



/(C/,a,/3) 



(27r)2 (27r)2 



M.2 ■ M.1 



[M - (Si - B^f] exp [^aBl - l3{Bi - R^f] 



It is evaluated in Appendix A with the explicit result (AS) 
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16n a^{a + P) 



(13) 



(14) 



From the induced momentum spectrum (11), (|T^) one may rederive the cj-spectrum per unit length ujdl /dujdz for 
soft gluon emission off quark jets produced inside the medium. After performing the exponential [/-integration and 
the r, Ti integrals explicitly the result of III is reproduced, namely 



ttL 



diodz 



In I cos WoTqI 



(15) 



We note from (||) and ( p^ ) that luqTq is expressed in terms of medium-dependent quantities and the gluon energy 
by 



UJ 



I qL /u!, 



2Cf 



(16) 
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where we introduce the transport coefficient 



q^^^P I dQ'Q'^. (17) 



Integrating (|T^) over u) the energy loss per unit length is obtained as 

- dE/dz ^ r ^du^^ qL, (18) 



so long as i? > Ecr — qL 
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III. INDUCED RADIATIVE ENERGY LOSS OF A HARD QUARK JET IN A FINITE CONE 



Let us consider a typical calorimetric measurement. We have in mind a hard quark jet of high energy E produced 
by a hard scattering in a dense QCD medium and propagating through it over a distance L. The quark loses energy 
by gluon radiation induced by multiple scattering. In the following we calculate the integrated loss outside an angular 
cone of opening angle 6'cone (Fig- 3), 



AE{9, 



duj 



jdl 



duodzdO 



■de. 



(19) 




gluon 



E ^ 


^^^^^l,-"""'^^ cone ^ 




quark 






E-AE (e,„„e ) 



FIG. 3. Hard process producing a quark jet. The gluon is emitted outside the cone with angle 
We note that for ^cono = the total loss ([l8| ) is obtained, namely j|] 



AE 



-qL . 



In detail we consider the normalized loss by defining the ratio 

Ai?(6'cone) 



-^(^conc) 



l\E 



which may be decomposed into 



cone ) — ^on- shell (^'cono) + R 



vertex l,"conc ) 



according to the treatment in the previous section. 
Defining C/cone i'^ the approximation of small cone angles 



■cone 2 cone' 



the B integral in ([l9| ) may be performed using 

f°° 1 
Jul ina^ 



(20) 



(21) 



(22) 



(23) 



(24) 
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The remaining u, r and ti integrations cannot be done analytically. It is convenient to change variables by introducing 
dimensionless ones: 



T = Toy, 

Tl = TqZ, (25) 
ujqt = (1 + i)x, 

which implies that the gluon energy is expressed by 

u;=^^\t,'vr^iy/x)'. (26) 

As a nice consequence, the ratio i?(0cone) turns out to depend on one single dimensionless variable 

R R{c{L)9 

cone ), (27) 

where 

c^{L) = ^q{L/2f. (28) 

The "scaling behaviour" of R means that the medium and size dependence is universally contained in the function 
c(L), which is a function of the transport coefficient q (^^ and of the length L, as defined by (p8|). For fixed L the 
medium properties are exclusively described by q which is different for cold and hot matter (cf. our discussion in 
Sect. 5 of II). The fact that ^cone scales as \/c{L) may be understood from the following physical argument the 
radiative energy loss of a quark jet is dominated by gluons having uo ~ gL^ (cf. Eq.(6.9) in I). The angle that the 
emitted gluon makes with the quark is 9 = k/oj with k the gluon transverse momentum. But k^ ~ qL so that the 



typical gluon angle will be 



1 



(9) 



3/2 



1.0- 




gl . \ . 1 . 1 . 1 . 1 

20 40 60 80 100 

C (L)e 

cone 

FIG. 4. Fractional induced loss R{Ocona) (solid curve) as a function of c{L)9cona, and decomposed into _Ron-shcii (dashed 
curve) and -Rvertox (dotted curve). 

For completeness the explicit dependence of R on c(L) ^conc is given in Appendix B. The numerical evaluation of this 
dependence is shown in Fig. 4 The following remarks may be useful. Obviously by definition i?(6'cone = 0) = 1, 
where we know from III that 

-Ron-shell (^'cono = 0) = 1/3 and -Rvcitcx(^'cono = 0) = 2/3. (29) 

At large angles, c{L)9conc ^ 1, we find (see Appendix C) the power behaviour 
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mf) ) ^ 4r(i/4) 1 

coming exclusively from _Ron-sheU- ^vertex (^cono) vanishes faster than the dependence given in (|30|), as may be seen 



from (C6) 



The ratio R{Oconc) is also universal in the sense that it is the same for an energetic quark as well as for a gluon 
jet, since the function c{L) depends on the product Ci?A = NcXgiuoni = CnXn for colour representation R). It is the 
integrated loss, however, which is bigger by the factor Nc/Cp for the gluon than for the quark jet. 

The "narrowness" of the angular distribution can be read off from Fig. 4: e.g. R{Oconc) decreases from 1 to 0.4, 
when c{L)9cone is increased from zero to 10. Explicit values for c{L) are discussed in Sect. 4 when comparing hot and 
cold QCD matter for fixed L 

c(i) _»c(L)|^Oj^j3, (31) 
when considered as a function of ^conc- 



so that R(0 cone) is much narrower than R(0 cone) 

HOT 



COLD 



IV. PHENOMENOLOGY 

We know from I - III that the transport coefhcient q ( |l7| ) controls the medium induced energy loss (|2^) as well as 
the p_L-broadening of high energy partons: both quantities are directly proportional to q. The medium dependence of 
the distribution i?(0conc) in ( |27| ) (see also Appendix B) is determined by the coefficient c{L) which according to ( [2^ ) 
depends on the square root of q. In this section we estimate the magnitude of c(L), with the aim of finding differences 
in the angular distributions R{Oconc) for cold and hot matter. In the following we will confirm the statement given in 
( ^ ) quantitatively. In order to be able to compare with our discussion already given in II (and in III) we keep the 
same numerical values for the quantity q as used there. 

(i) For hot matter at reference temperature T — 250 MeV we take 

g~0.1GeV^ (32) 



This value is deduced from (17) using the screened one-gluon exchange cross section, e.g. for two flavours one derives 
§ 

AO 

q~-ai)valT^ ({3)^ 1.202. (33) 

TT 

We use a typical value for the dimensionless quantity w ~ 2. Eq.(|3^) shows explicitly the temperature dependence of 
q. For the values of T under consideration (^2|) leads to a rather large value of c{L), namely 

c(L) ~ 40(L/10/m)^/% (34) 

HOT 

where we use — 1/3 as in II. 

This large value of c{L) implies that for cone angles 6*00110 > 30° the normalized distribution R{9conc) is well 
approximated by its asymptotic form given by (pO|). For large c{L)9conc the fraction of the induced energy loss 
R{c{L)9conc) decreases as (TL)^'^/'* for increasing temperature T and length L (for fixed 6'cono)- We also note that the 
energy loss outside the cone quickly drops with increasing Scone (cf. Figs. 4 and 5): for ^conc — 10° (40°) the fraction 
R is reduced to 40% (20%) of the total loss, for which we find a value of 



AS ~ 60 GeV {L/\Qfmf (35) 



for a quark jet, according to Eq. (po|) . 
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0.2- 



10 20 30 40 

9 (degrees) 

cone ^ 

FIG. 5. Medium induced (normalized) energy loss distribution as a function of cone angle 6'cone for hot (T = 250 MeV) (solid 
curve) and cold matter (dashed curve) at fixed length L — 10/m. 

Indeed the medium dependent distribution R{9conc) is rather narrow on the partonic level (Fig. 5). 
(ii) For cold matter we quote from II the typical value 

g~ ^a,[a;G(x)]GeVV/m 
25 

~ 0.005 GeV^ 

with a gluon distribution, xG{x) ~ 1 — 2. Compared with ( ^5|) we find 



A£:~ 3.5GeV {L/lOfmf 



and accordingly a much smaller value for 



ciL)\ 



COLD 



10 {L/lOfmf^^ . 



This confirms Eq.(pl^) quantitatively, namely 



ciL) 



HOT 



4c(i) 



COLD 



(36) 



(37) 



(38) 



(39) 



at T = 250 MeV. When comparing the energy loss in hot and cold matter we observe from (|35| ) and (^^ that the 
induced loss in cold matter is much smaller (by a factor of about 15) than in hot matter, and that i?(^conc) 



IS 



broader than i?(0conc) 



COLD 



HOT 



which can be seen from Fig. 5. 



As a conclusion we expect that the medium-induced AE(6conc) for energetic jets can be large in hot QCD matter 
and although collimated, still appreciably larger than in cold matter even for cone sizes of order ^conc — 30°. 

So far we have discussed the medium-induced energy loss and its angular distribution. Concerning the total energy 
loss of a jet of a given cone size it is important to take into account the medium independent part. It corresponds to 
the K = contribution, which we have subtracted in Eq.(H). For the case of a fast quark produced in the medium the 
a = soft gluon radiation spectrum is given in leading order by the well known Born term. 



Ujdr=° UsCf 1 



(40) 



valid for single gluon emission off a quark in the vacuum, and therefore truly medium independent. This result may 
be obtained as follows: the k = factorization term can be realized as the emission of large oj gluons. These gluons 
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are produced with large formation times, such that the dominant gluon emission takes place after the energetic quark 
has propagated through the medium, namely far outside the medium. A more careful analysis confirms the result 
(^). The angular dependence of the corresponding energy loss of a quark jet in a cone may be estimated from (|4C|), 



Lodr 



K=0 



diodzcPU 



asCp „ , / ^max 

2 A in 



(41) 



using a constant Og. Since in this case the loss is sensitive to gluon energies uj = 0{E), one should take into account 

the quark gluon splitting function f ^'''^ V^"* ' which amounts to replacing E by 2/3E in 

Keeping in mind that the following estimates are based on the leading logarithm approximation, we summarize in 
Fig. 6 the medium-induced (for a hot medium with T = 250 MeV) and the medium independent energy losses. 



> 
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FIG. 6. Energy loss in a hot medium, T = 250MeV, as a function of ^conc- Dashed curve represents ( |4l| ) for E = 250GeV. 
L = 10/m. 

For the strict validity of we recall the condition Q that E > i?crit — qL'^ ^ which requires very large values 
of the quark energy: E > 250 (L/lO/m) GeV, when T ~ 250 MeV. Such energetic jets may indeed be produced in 
hard scatterings in heavy-ion collisions at the LHC at CERN. For L — 10 fm this jet energy leads to Aii^''^°(6'cono) — 

15 In g"*" ^ GeV. To get this estimate we have taken ag — 1/10 due to the large characteristic transverse momentum 

of the emitted gluon. The ^cone dependence is shown in Fig. 6 for the two contributions: over the limited range of 
10° < 6*00110 < 40° they turn out to be not very different (^max is taken C'(7r/2)). It is, however, important to keep in 
mind that ( ^ ) and ( |20| ) differ with respect to their dependence on L. From ( |4l] ) AE'^^'^{0cone) does not depend on 
L at all. The L dependence of Ai?(6'conc) = ^(6'conc)Ai? comes from two sources. It is proportional to from AE, 
with an extra L dependence contained in R{9conc) at fixed cone angle due to c{L), Eq.([28|). 
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APPENDIX A: B-SPACE INTEGRATION 



The gluon momentum spectrum is expressed in terms of an integral I{U_, a, f3) (Eq.([l^)). It can be reduced to 



I{U,a,P) 



d_ d_ 

d[3 da 



J{U,a,P), 



where 



J(C/,a,/3) 



(fBi d?B2 ,(B,- 



(27r)2 (27r) 

^2 



X Si • =1 exp [-aB? - I3{B^ - B^f] . 



B.2 



In order to perform the B_i , B_2 integrations we note that 



B2 I 
bI~2 



Vb, In Bl, 



and 



B, exp [-aBl - PiB, - B^f] = -— ( Vs, + Vs, ) exp [-aBj - (3{Bi " R^)^] ■ 



These identities allow us to write 



JiU a 0)^~— I ^^-M^e'iB,-B,yu 



Vs. In Bi^ (Vs, + Vs.) exp [-aBj - P{B^ ~ B^ f] 



After partial integration and using 



V|, InBl^^TrS^B^), 



one finds 



J{U,a,f3) = 



1 



1 



167r2 a(a + /3) 



exp 



4(a + 13) 



Finally the function I{U_, a, (3) becomes 

J(C/,a,/3) 



1 



1 



167r2 a'^{a + l3) 



exp 



4(a + /5) 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



APPENDIX B: EXPLICIT EXPRESSION FOR i?(6lcoNE) 



Here we summarize the explicit expressions for i?on-shcii(^conc) and i?vortex(^conc}- They are obtained from Eq.(|l] 
performing the integrations defined in (^9|) by applying the change of variables (| 
The results are: 



Ro 



sholU 



4 

Re 



dx 



dy 



dz < exp 



2 

cone 



(l-z) 



((l + ^)^f 
sin^(l + i)x 



exp 



2 

cone 



tan(l+^)a; 
y~(TTi)x~ 



(Bl) 
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and (after shifting z^\/y~\~z /y) 



R- 



vertex I, L'conc ) 



dx 



ydy 



dz 



(z + iy 



cxp 



2 

cone 



Re 



cos2(l + i)x 



z + 



taii(l+'i)a: 
{l+i}x 



exp 



(l + yz) \x 



92 

cone 



y (iW +^ + y^-y 



(f)* 



(B2) 



The subtractions in ( |B1| ) and (|B2| ), the k — term in (g), are taken at fixed U_ rather than at fixed 0^^^^, as is 



necessary for the subtracted term to be medium independent. 



APPENDIX C: BEHAVIOUR FOR LARGE 6Icone 



From the expressions (Bl) and (B2) one may evaluate the behaviour of R{Oconc) at large values of c(L)0cone ^ 1- 
In view of the exponential dependence on c{L)9conc of the integrands we conclude that this behaviour is controlled 
by the endpoint in the variable y, y = 0. Expanding the integrand in (Bl) arround this point we find 



4 f°° dx 

-Ron- shell (^'conc) ^ —Re 



1 + i)x 
sin(l + i)x 



X I dy dzexp 
'0 Jo 



(l-z) 



eone / y_ 
X 



Performing first the Gaussian y^-integration and then the z-integration we obtain the asymptotic behaviour 

2' 



R 



•on — shell l,C7cone 



1 



4 r(i/4) 

4 r(i/4) 1 

5 ^ 7c(^)^ 



1 



(1 + 
sin(l + i)x 



The x-integration is explicitly given by noting that 



Re 







dx — 

dx 



(1 + 



1 



x tan(l + i)x 



The corresponding expansion around j/ ~ in ( p32D leads to 

r-l 



R 



■vertex I, "eone 



c(L)eoo„o»l TT 
1 



4 dx 
-Re 



ydy 



dz 



{z + iy 



cos2(l + i)x 



z + 



tan(l+i)a; 



exp 



-C^WCne(f)' 



where the z-integration can be performed immediately. The y-integral leads to 



dy^ exp 



4 c{L)ec 



■ erf 



where we keep the erf-function ||l8| in order to have a finite x-integral at fixed, but large c:{L)9c 



Rv 



c(^cone) 



1 



1 



c(L)9,onc»l c{L)d, 

X Re 



X 



X^ 



{l + i)x 



sin(l + i)x cos(l + i)x 



(CI) 



(C2) 



(C3) 



(C4) 



(C5) 



(C6) 
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In Fig. 7 we show that the asymptotic behaviour of (C2) and ( |C^ ) is also confirmed numerically [Q. The large 
^conc~behaviour of i?(0conc) is dominated by i?on— shell (^conc 

) given by (|C 
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FIG. 7. -R nn-s hrii (dashed-dotted) and i?vcrtGx(Sconc) (solid) compared with the asymptotic forms (dashed and dotted curves) 
as given by (C2) and (Ct), respectively. 
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